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Abstract 

We consider integrable models of the Haldane-Shastry type with open boundary con- 
ditions. We define monodromy matrices, obeying the reflection equation, which generate 
the symmetries of these models. Using a map to the Calogero-Sutherland Hamiltonian of 
BC type, we derive the spectrum and the highest weight eigenstates. 
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It has recently been understood that while the long-range interacting spin chains 
capture some of the essential physical features of the Heisenberg spin chain, exact results 
can often be obtained more explicitly. In particular, the spectrum can be described in 
terms of elementary excitations called spinons. The expressions of the wave functions can 
be derived. In this paper we propose generalizations of the Haldane-Shastry Hamiltonian 
@ which describe open spin chains with various kinds of boundary conditions. 



1. The Hamiltonians 

We consider a system of spins on a one-dimensional semi-circular lattice, with an 
interaction depending on the inverse square of the distance between the spins 

N N 

H = Yl N(^j-l)+^(^-l)]+^(&i^o + 26 2 ^)(^-l). (1.1) 

v£j=\ i=l 

The operators P\j permute the spins at the sites i and j. They correspond to a Heisenberg- 
type interaction, P^ = (<rfcr° + l)/2, with a a the Pauli matrices. Pi is a reflection operator, 
Pf = 1, whose representation will be chosen later and P^ = PiPjPij. The coupling 



constants are defined as follows: 



ZiZj - _ ZiZj 
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hij — , * 2 , hij — hio — / 1N2 (1-2) 

(zi-Zj) 2 (Zi-Zj ) 2 (zi-1) 2 

with Zi points on the unit circle characterizing the positions of the spins. 

There are at least two possible representations for the reflection operators Pf, the 



simplest one is to put Pi = 1, which leads to Pij = Pij. The Hamiltonian ( |1 . 1|) then 
becomes 

N 

H= (hij + hij) (Pij -1). (1.3) 

This model is sit (2) symmetric. It was first considered by Simons and Altshuler 0. 

The reflection operators can be equally represented by the Pauli matrices Pi = cr|, in 
which case the sit(2) invariance is lost; only the projection of the total spin is conserved. 

These models can be interpreted as open versions of the Haldane-Shastry spin model 
0. The terms involving hij couple the 'real' spins % and j, while the terms with coupling 
constants hij correspond to the interaction of the spin % with the image of the spin j. 

When Pi = erf, the last two terms in (|1.1|) can be interpreted as space dependent 
magnetic fields varying as l/sin 2 x, with x being the angle (or two times the angle) between 
a spin and the boundary. Remark that for the term ha the two boundaries are equivalent, 
while for the term hio they are not. This term breaks the z% <-» — z% invariance. 
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2. Integrability Condition 

This Hamiltonian is integrable only for certain values of the numbers z% and of the 
parameters 61, 62- To obtain the condition of integrability we shall follow the approach of 
ref. HJ]. In addition to the operators which permute the spins , let us introduce permutation 
and reflection operators acting on the coordinates as 

Introduce also a set of mutually commuting operators (Dunkl operators) || 

d t =J2 9 H K ij ~ 6 3 iKi + Y^^^ + + b ^n) K i ( 2 - 2 ) 

j;j>i j;j<i i^j 

with 6ij = Zi/(zi — Zj) : Oij = Zi/{zi — zj 1 ) and 6io = Zi/(zi — 1). We denote by n a 
projection operation which consists in replacing the coordinate permutations (reflections) 
by spin permutations (reflections) when they are at the right of an expression. Let us 
consider the Hamiltonian H with the property H = ir(H) 

N N 

H= Yl [M^-^ + M^i?- 1 )] + ^(b 1 h i Q + 2b 2 h ii )(K i -l). (2.3) 

ii=3=X i=l 

This Hamiltonian commutes with the reflection and permutation operators, [H, Kj\ = 0, 
[H, Kij] = 0, so its eigenfunctions can be chosen either symmetric or antisymmetric under 
the reflection z% — > z~ x and permutation of coordinates. H can be diagonalized on a basis 
of functions depending on coordinates and spins, These are also the eigenfunctions of 
H, provided that (K t - P-i)^ = 0, (K^ - P^)* = 0. Choosing K { = +1 or -1 lead to the 
choice Pi = 1 or Pi = af. 

The integrability condition is that the Hamiltonian (|2.3| ) commutes with the Dunkl 
operators, [H : di] = 0. It can be translated into iV equations for the coordinates z$ and 
the constants b\,2 

N 

2 Y ( h ij w ij + h ij w ij) + 2 hiiWa + 61(61 + b 2 ) h iQ Wio = (2.4) 

where Wij = {zt + Zj)/(zi - Zj), Wij = {zt + z~ 1 )/(z i - zj 1 ) and w i0 = (z t + l)/(zi - 1). 
The system of equations Yl^-j^i^i^ij = nas the solution Zk = Ce 2l7rk / L with C a 
multiplicative constant. We derive solutions for the equations ( |2.4|) by putting them in 
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i) L=2N ii) L=2N+2 L=2N+1 

Fig. 1: Distribution of the spins on the lattice. The bold line represents the 
boundary (mirror), the filled circles indicate the positions of the spins. The 
open dots are the images of the 'physical' points through the mirror. 

this form. Different solutions correspond to different distributions of the coordinates Z{, 
as illustrated in Fig. 1. 

We consider a circular lattice with L sites and set uj = e 2in / L ; 

i) If the 'mirror' (the boundary) is placed between the sites of the lattice, then L = 2N, 
N being the number of 'physical' sites. In this case, Zk = uj k ~ 1 / 2 and the constants in 
( Pn|) have the values bi = 0, 62 = ±1- 

ii) If the mirror passes through the sites of the lattice, L = 2N + 2 and Z\~ = u k . The 
constants in the Hamiltonian have the values b± = 0, 62 = ±3. 

Hi) The third case corresponds to a mirror passing through a site and a bond, so L 
is odd, L = 2N + 1. The positions of the sites are Zk = uJ k and 62 = ±1- The constant b\ 
acquires a non-zero value and it is the solution of the equation 61(61 + 62) = 2. 

When the reflection operator is represented by the identity, the constants b\ and 62 are 
not relevant, as they multiply constant terms in (|1 . 1| ) . When the reflection is represented 
by of, an overall change of sign of the constants b±, 62 corresponds to the conjugation 
of — > — of, of — > of. In this case we will retain only the positive value of 62- 

3. Symmetry and Spectrum 

The model possesses a monodromy matrix which commutes with the Hamiltonian H 
and which satisfies the reflection equation || : 

R w {u - v)T°(u)R 00 ,(u + v)T°\v) = T '(v)R Q o'(u + v)T° (u) R 00 > (u - v) (3.1) 
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where T°(u) stands for T(u) <S> 1 and T° (u) for 1 ®T(w). The matrices R(u) and R(u) are 
given by : 



R(u) =u + Pi 



00' 



R(u) = u + P 



00' 



(3.2) 



where the operator P o' permutes the two auxiliary spaces and 0' and Poo' represents per- 
mutation followed by a reflection in each space. The explicit expression of this monodromy 
matrix is 



T(u) 



TV 



N 



n i+ 



Pi 



u — di 



1 + 



h + b 2 P 



u 



--N 



n i+ 



p 



t0 



u + di 



(3.3) 



Expanding T(u) in powers of u 1 around the point at the infinity generates conserved 
charges, T(u) = I + X] n >o T n u~( n+1 \ If Pi = 1, the first conserved charge is the total spin 



Qq = YliLi a t anc ^ ^ ne non-trivial charges are located at even levels, Q2n- If Pi = of j the 



conserved charge at the level zero is the projection of the total spin, Qq — ^2^ =1 erf. The 
su(2) invariance is broken by the presence of the last two terms in the Hamiltonian ( |1 . 1| ) . 
At the next level, we find two non-trivial generators : 



Qf 



i=i 



^ {wijPij + WijPij) + hwioPi + b 2 waPi 



(3.4) 



As a consequence of these symmetries, the energy levels are degenerate and the eigen- 
vectors are grouped in multiplets. Each multiplet is characterized by a set of rapidities 
{mj, C1/2 + C2 < nii < L/2 — 1, where ci, C2 and L/2 — rrii are integers. The values of c\ 
and C2 corresponding to the different cases are summarized below 







Pi 
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Pi 








L 


rrii 


Cl 
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Cl 
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Table 1. 

The energy of a multiplet have the expression 



E 



{mi} 



M 

E 

3=1 



T 



(3.5) 



The form of the spectrum will be justified later, using a map to a continuous model, the 
Calogero-Sutherland model. 

The rules giving the degeneracies of the spectrum depend on the symmetry of the 
model, hence on the representation chosen for the operators Pi. They have been determined 
numerically and they should also follow from a careful analysis of the representation of 
the symmetry algebra ( |3.1| ). If -Pi = 1, the model possesses a sit (2) symmetry and the 
degeneracies are the same as for the Haldane-Shastry model 0]. The only difference is 
that the energy is not invariant under the transformation m; — > N — rtii. The rule can 
be formulated as follows: take a chain of iV fictitious spins with values + and — (these 
are not the spins in (|1 . 1|) ) . Put a 1 between two consecutive fictitious spins if the first 
one is greater than the second and a otherwise. This sequence of and 1 is called a 
'motif and the positions on the motif are labeled by integers (half-integers) ranging from 
ci/2 + C2 to L/2 — 1. The positions of the l's give the set of quantum numbers {mi}. The 
degeneracy associated to a given motif equals the number of configurations of the fictitious 
spins compatible with it. 

If Pi = of, the symmetry group is w(l) and we expect the spectrum to be less de- 
generate. In this case, the rule giving the degeneracies is slightly modified from the one 
above. We must consider now a chain of N + 2 fictitious spins, with the first and the last 
ones fixed at the value +. This constraint imposes to all the fictitious spins at the left 
of the first 1 to have the value +. The positions of the symbols on the motif range from 
ci/2 + C2 to L/2. For the case L = 2N + 1, b\ = —2 we obtained numerical evidence for a 
su(2) symmetry, but we did not succeed to find the conserved generators (which are not 
the total spin operators). In this case, the degeneracies are the same as for the periodic 
chain. 

4. Mapping to the Calogero-Sutherand Model 

The eigenvectors and the eigenvalues of the Hamiltonian (|1.1|) can be obtained from a 
correspondence with a particular model of Calogero-Sutherland (CS) type. Before giving 
the details of the correspondence, let us give the definition of the general CS model. This 
model describes n particles on a line, with coordinates Xi, < x\ < I. In the version 
proposed by Calogero and Sutherland M || , the interaction among the particles is pairwise 
and the model is completely integrable. In fact, the model is integrable for more general 
potentials, associated to all the root systems of the simple Lie algebras [[Kj . 
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Let V denote a n dimensional vector space with an orthonormal basis {ei, e n } and 
let R = {a} be a root system spanning V, with R+ the set of positive roots. Let x denote 
the vector (xi, x n ) and x ■ a its scalar product with the vector a. Then the general CS 
Hamiltonian is 

where g a is constant on each orbit of the Weyl group (it has the same value for the roots 
of the same length). The Hamiltonian proposed by Sutherland and used to solve the 
Haldane-Shastry model correspond to the root system of type A n , R = {ei — ej, i 7^ j}. 
The one relevant for our problem is associated to the BC n (generalized) root system, 
R = {±ej, ±2ej, ±e; ± ej}. In this case, the Weyl group is generated by the n reflections 
on the hyperplanes Xi = (Ki) and by the reflections on the hyperplanes Xi = ±Xj (the 
permutations Ky, Kij). 

With the variables Zk = e lXh the Hamiltonian (|4.1|) becomes: 

n n n 

H C s = ^(^a 2 J 2 +/3(/3-l) ( h H +E M 2c2 + Cl - l ^ hi0 + 4c2 ( C2 " ■ 

j=l i^j=l i=l 

(4.2) 

with hij defined in terms of z as in ( |1 . 2j ) . The form of the three independent coupling 
constants g a was chosen for later convenience. The ground state wave function of this 
Hamiltonian is 

Mz) = n^" (/3(7t " 1)+C2+Cl/2) (^ - l) Cl k 2 - 1) C2 \\izi - ztfiziZj - If (4.3) 

i i<j 

In some situations, a gauge transformed version of (|4.2| ), Hcs = <Po 1 Hcs4>0i can be useful 

n n 

Hcs = ^2(zid Zi ) 2 + (3 ^(wij + Wij)zid Zi + ^(ciu> i0 + 2c 2 wn)zid Zi + E (4.4) 

i=l i^j i=l 

with Eq = Y^i-i [ c i/2 + C2 + (3(n — i)] 2 . As in the case initially studied by Sutherland, 
this Hamiltonian is triangular in a basis of symmetrized plane waves 



J> Al 4 2 -*n" (4-5) 



where the sum is now over all the permutations and the changes of sign of the integers Xi 
(Ai > ... > A n > 0). 
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Therefore, the eigenfunctions of the Hamiltonian Hcs can be taken of the form 



4>(z) = Uz)PM,...x n {z) (4.6) 

where the P(z) are polynomials in the variables Z{, z~ , symmetric under the permutations 
Zi <-> Zj and under the reflections Z{ — > 2~ . Some properties of these polynomials associ- 
ated to general root systems are studied by Macdonald jnj . The corresponding eigenvalues 
are 

n 

Ecs = ^(A, + 0(n -i) + a/2 + c 2 ) 2 (4.7) 
i=i 

The spectrum and the eigenvectors of the spin chain can be derived from the spectrum 
and the wave functions of the CS model . We construct the Hilbert space of the chain 
from the ferromagnetic state \Q) = | + + ...+) by reversing M spins : 

N 



l*)= ^,...,n M (T ni ...CT nM \n) (4.* 



ni,...,n M = l 

The coefficients ipm,...,n M are symmetric under the permutation of their indices. We extend 
their definition for 1 < n < L by putting V'...,L-n i +2a,... = Ki ifi,.. jni ,... = ±V'...,ni,... with 
a = 1/2 if L = 2N and a = otherwise. The + sign corresponds to the choice Pj = 1 and 
the — sign to Pi = af. By convention, we chose the coefficients equal to zero if an index 
equals L or if two indices coincide. 

Let now Zm) be a symmetric polynomial of degree at most L — 1 in each Zi 

such that i>m,...,nM = ^(uj ni ~ ct , ...,uj nM ~ OL ). The condition of antisymmetry of ip (when 
Pi = af) and the condition for to be a su{2) highest weight, Ym=i a t\^) = (when 
Pi = 1) impose that =1 ^ ni ,...,n M = 0- This means that ^{z\ = 0, Zm) = and 

the degree of *&(z) in each variable is greater than 0. 

We want to translate the action of the spin Hamiltonian on the vectors \^>) into an 
action on the functions ^f(z). To do that, we use the fact that in the basis of polynomials 
Qk in one variable of degree (L — 1) specified by Qk{oj n ~ a ) = 5%, the matrix elements of 
the derivatives are: 

- ^)Qk{z) = - E -J^Qii*) 



(4.9) 

zd z (zd z - L)Q k (z) = -2 J] , k^'j^ (Qj( z ) - 



L U) k U) j 



This enables to recover the action of the Calogero-Sutherland Hamiltonian on the functions 
H\V)= (-!)^ ni 

ni ,...,n M 

The values of the coupling constants in Hqs are the values given in Table 1 and the 
constant (3 is equal to 2. Note that, for these values, the energy in fl4.7p is consistent with 
(pT5|) . The rapidities m; are defined in terms of the integers Aj as m; = \ + 2(M — i) + 
ci/2 + C2- They are integers or half-integers, depending on the parity of ci, and they obey 
the following selection rules : m 8 — rrii+i > 2 and c±/2 + C2 < rrii < L/2 — 1. The bounds 
on the rrii follow from the fact that ^(z) is a polynomial of degree between 1 and L — 1 in 
each variable z%. 

Thus, we have obtained a set of eigenfunctions of the spin chain from the wave func- 
tions Q4.6Q of the Calogero-Sutherland model. This procedure gives only the highest weight 
vector (annihilated by Qq or Q^) in each multiplet. The other vectors in the multiplet 
can be obtained by repeated action of the generators Qo,i- 
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Hcs- 



,00 



n M —OL 



)a-,...,a- M |0) (4.10) 
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